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We set up a structural model to study credit risk for a portfolio containing several or many 
credit contracts. The model is based on a jump-diffusion process for the risk factors, i.e. for the 
company assets. We also include correlations between the companies. We discuss that models of 
this type have much in common with other problems in statistical physics and in the theory of 
complex systems. We study a simplified version of our model analytically. Furthermore, we perform 
extensive numerical simulations for the full model. The observables are the loss distribution of the 
credit portfolio, its moments and other quantities derived thereof. We compile detailed information 
about the parameter dependence of these observables. In the course of setting up and analyzing our 
model, we also give a review of credit risk modeling for a physics audience. 
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I. INTRODUCTION 

Economics attracts the interest of a quickly growing 
community in physics. A large part of the research ad- 
dresses the financial markets. Attempts are being made 
to better understand various phenomena such as the fat 
tails of the stock price distributions by relating them to 
physics systems, see Refs. [TJH1I2] and references therein. 
Physicists have also joined the activities of economists 
and computer scientists in agent -based models [U [5], 
and applied their long-standing experience in complex 
systems, see Ref. [BJ. 

As far as the field of finance is concerned, the vast 
majority of studies put forward by physicists has been 
devoted to market risk. The market risk is due to the 
unknown time evolution of the asset prices. In general, 
one is faced with a large spectrum of different risk types. 
One also distinguishes, for example, operational risk (due 
to failure of internal systems), political risk (due to po- 
litical decisions that affect the capital markets) and legal 
risk (due to fraud and discontinued contracts). In this 
contribution, we address credit risk. It is due to the fail- 
ure of a counterpart to make a promised payment. At 
present, risk managers and researchers are more familiar 
with market risk than with any of the other risk types 
and the corresponding mathematical description is highly 
developed. It is of considerable practical interest to im- 
prove the knowledge about and the modeling of the other 
risk types. In the case of credit risk, the probability that 
a promised payment is not made is usually small and 
difficult to estimate. Nevertheless the amount of money 
involved and thus the associated loss can be enormous 
and even jeopardize the existence of the financial insti- 
tution which issued the credit. 

Only recently, physicists started applying their spe- 
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cific tools to credit risk HI E EU EEU E2]- The inter- 
esting point for practitioners and researchers, especially 
statistical physicists, is the highly asymmetric form of 
the loss distribution and the resulting peculiar features. 
This distinguishes credit risk from market risk, although 
the former clearly depends on the latter. In investment 
theory the standard deviation, referred to as volatility, 
of the relative asset price change is taken as a measure 
of how risky a certain investment is. If more uncertainty 
is incorporated in the investment, i.e. if the volatility is 
larger, then the demanded earnings are higher. Due to 
this fact, investors are traditionally risk averse. In other 
words, a potential loss is considered to be more puni- 
tive than a potential gain is beneficial, even if they are 
equally probable and large. The asymmetric character of 
the loss distributions makes risk measures other than the 
volatility also important in credit risk management. 

In this study, we set up and analyze a structural model 
for credit risk, based on a jump-diffusion process for the 
risk factors. Our study is related to, but different from 
the work in Ref. 1 13] . These models are particularly ap- 
pealing to physicists, because their starting point is, in 
physics terminology, microscopic and dynamical. This 
gives them a rather general character which makes them 
also suited for other problems in physics and in the theory 
of complex systems. With this contribution, we pursue 
two goals. First, we systematically explore the interplay 
between the different parameters of our structural model, 
particularly the role of leverage, jumps and correlations. 
In contrast to the existing literature, our main focus is 
on the full loss distribution of the credit portfolio, its mo- 
ments and its tail behavior. Second, we review credit risk 
modeling and keep the whole presentation pedagogical, 
because we want to make this topic more accessible to 
the physics audience. 

The paper is organized as follows. We review the 
present status of credit risk modeling in Sec.[n] In Sec. Ill 
we introduce our model. We discuss a simplified version 
of it analytically in Sec. |IV| and the full model numer- 
ically in Sec. IV] Summary and conclusion are given in 



2 



t=0 



t=T 



Creditor 



Principal 



Obligor 



Creditor 



Obligor 



FIG. 1: The cash flow of a zero-coupon bond. 



there has been a missed or delayed payment of interest 
and/or principal within a grace period, or that an obligor 
files for bankruptcy [16 . Although default is a truly rare 
event, creditors can lose large amounts of money. A good 
example is the bankruptcy of Enron Corporation in Dec. 
2, 2001. At an asset value of US $49.53 billion, this 
was the largest bankruptcy filing in US history to that 
date [IT]. The actual loss for the creditors was US $9.9 
billion [IS]. 



Sec.|VH 



C. Credit ratings 



II. CREDIT RISK MODELING 



After defining debt instruments in Sec. |II A| we discuss 
what one means by default and credit ratings in Sees. II B 
and |II C| respectively. The credit risk measures are intro- 
duced in Sec. |IID| The impact of correlations is discussed 
in Sec. |IIE| The role of the probability density function 
for credit losses is explained in Sec. II F We conclude this 
short review by presenting the most important credit risk 
models in Sec. Ill Gl Reviews in the financial literature on 
credit risk modeling can be found in Refs. [T3J H3] 



A. Debt instruments 

A debt instrument is simply a written promise to re- 
pay a debt. There is a wide range of such contracts. A 
debt instrument has two positions: a lending side (the 
creditor) and a borrowing side (the obligor). Bonds are 
very common. A bond is issued for a period of one year 
or more with the purpose of raising capital by borrow- 
ing. The government, states, cities, corporations, and 
many other types of institutions can sell bonds. Gen- 
erally, a bond is a promise to repay the principal along 
with interest (coupons) on a specified date (maturity). 
The principal is the amount borrowed or the part of the 
amount borrowed which remains unpaid (excluding in- 
terest). Some bonds do not pay interest. In this study 
we focus on the zero-coupon bond. The cash flow of the 
zero-coupon bond is limited to two dates: the date of 
issue t = and maturity t — T. At the issue date the 
creditor lends a specified amount of money to the obligor. 
At maturity, the obligor has to repay the face value of the 
bond. The face value is the amount borrowed plus inter- 
est and additional yield compensating for the risk. Fig.[l] 
shows the cash flow of the zero-coupon bond. 



B. Default event 

The key issue which separates credit risk from, for ex- 
ample, market risk, is the concept of default. Its def- 
inition is not homogeneous throughout the industry of 
credit risk management. Usually default means that 



A way to quantify credit risk is to determine the credit 
worthiness of a potential customer from the historical 
performance of the obligors. There is a wide range of rat- 
ing systems for credit worthiness, systems used internally 
by the credit institute as well as external ratings available 
for public. The credit rating of a company is often di- 
rectly linked to the probability that the company defaults 
within a fixed time horizon, usually one year. Two fre- 
quently used rating systems are Standard & Poor's and 
Moody's. The probability that a company changes its 
credit rating is expressed in terms of a rating transition 
matrix which contains the probabilities that a company 
with a certain rating migrates to another category, usu- 
ally within a year. The credit rating transition matrix 
is based on the historical migration frequencies of corpo- 
rate bonds. It is observed that the most probable future 
event is that the company remains in the same rating 
category This is valid for all rating categories. More- 
over, the probability of a downgrade is generally higher 
than the probability of an upgrade [15] . 



D. Credit risk observables 

There are several ways of quantifying credit risk. 
We distinguish between standalone risk and portfolio 
risk |19j . The most frequently used standalone observ- 
ables are the default probability (DP), the loss given de- 
fault (LGD) and the migration risk. The conventional 
portfolio risk observables are default correlations and ex- 
posure, i.e. the size, or proportion, of a credit portfolio 
exposed to default risk. 



E. Correlations 

The credit worthiness of obligors often involves mutual 
correlations. For example, defaults are more frequent in 
times of regression in the surrounding economy. Further- 
more, one can see that companies in the same country 
and/or industry can affect each others rating migrations, 
up as well as down [14]. There are different ways to in- 
corporate correlations between default events of obligors. 
For example, one uses the correlations between the eq- 
uity values of the companies, i.e. the stocks, to describe 
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the dependency of credit migrations. Another way to ex- 
amine this is to look at how the obligors depend on the 
current state of the economy [5U] . Modeling correlations 
is often difficult because only limited data are available 
for the indicators chosen. Moreover, because of the huge 
number of correlations involved in a normal sized credit 
portfolio (which can contain a few thousand bonds), it 
is necessary to make simplified assumptions. A common 
procedure is to categorize the companies into different 
groups or branches with specified correlations. For ex- 
ample, one can assume correlations which are country- 
specific, industry-specific etc. Even though these simpli- 
fications lead to a more manageable model, it is still a 
complicated task to decide the structure of the catego- 
rization and to estimate the group-specific correlations. 



F. Probability density function of credit losses 

The primary output of a credit risk model is the prob- 
ability density function (PDF) of credit losses for a given 
portfolio. Adapting to the more common physics termi- 
nology, we refer to this function as the loss distribution. 
Here we make the assumption that the loss is a continu- 
ous random variable, so that we can work directly with 
the PDF. From the loss distribution one determines the 
expected loss EL, the unexpected loss UL, the required 
risk capital and further quantities. The expected loss EL 
is the mean of the loss distribution and the unexpected 
loss UL is its standard deviation. It is important to notice 
that UL, not EL, measures risk. However, to cover po- 
tential losses it is not enough to have a "cushion" within 
the standard deviation. The probability that losses ex- 
ceed the UL is significant and it is therefore necessary to 
have another measure of risk capital. To quantify risk 
capital one usually uses the economic capital EC which 
is also known as Capital at Risk CaR or as the Value at 
Risk VaR. The economic capital EC is defined as the dif- 
ference between the expected loss EL and the a-quantile 
for a certain level of confidence. Figure [2] shows a typ- 
ical loss distribution with EL, UL and EC. The general 
appearance of such a loss distribution is different from 
the distributions in, for example, market risk. While the 
distributions of market risk are typically Gaussian, the 
distributions generated by credit risk are skewed and lep- 
tokurtic [14] . 

To measure the important tail behavior, one uses the 
kurtosis excess. We recall that the kurtosis is defined 
as the fourth central moment divided by the squared sec- 
ond central moment. The kurtosis excess 72 is the part of 
the kurtosis that exceeds the kurtosis of the normal dis- 
tribution (which is equal to 3), i.e. 72 = fii — 3. A distri- 
bution is often referred to as fat -tailed if it is leptokurtic, 
i.e. 72 > 0. Fat-tailed distributions have higher quantilcs 
than a normal distribution and thus require more atten- 
tion from a credit risk manager. 

To generate a loss distribution of a credit portfolio, one 
can evaluate it for a structural model either analytically 
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FIG. 2: A schematic loss distribution with the expected loss 
EL, the unexpected loss UL and the economic capital EC. 



(in some simple cases) or numerically by Monte Carlo 
simulation. It is also possible to approximate an actual 
portfolio distribution by some known analytical distri- 
bution. In the latter case, one maps the actual portfo- 
lio with unknown distribution to an equivalent portfolio 
with a known distribution. A frequently used distribu- 
tion which has the right features, i.e. skewed and fat- 
tailed, is the beta distribution [14]. A better and more 
sophisticated approach than parametric approximations 
is offered by reduced form models. Some of these mod- 
els lead to analytically tractable loss distributions, see 
e.g. Ref. [2T] and references therein. 



G. Credit risk models 

Current credit risk models can be divided into two 
main categories: structural models and reduced form 
models. 

In the structural models one makes assumptions about 
the time evolution of the risk factors, i.e. mostly the asset 
or stock prices of the companies, as well as about the 
liabilities. Whenever the asset value falls below some 
specified threshold, like the book value of the liabilities, 
the firm defaults. The structural credit risk modeling 
approach has its roots in the Black and Scholes theory 
for option pricing, and the Merton model, see Ref. [T3] . 
The Black and Scholes theory is based on the assumption 
of a friction-less market where the stock or asset price 
is described by a geometric Brownian motion. Merton 
viewed the equity value of a company as a call option of 
the asset value with a strike price equal to the face value 
F of the debts. It is assumed that the company has 
a certain amount of zero-coupon debt due at a future 
time T. In consequence, it is possible to apply the whole 
Black and Scholes machinery to the credit risk modeling 
problem [T3]. This is the "microscopic" viewpoint which 
makes structural models suitable for physics approaches. 

The reduced form models for credit risk are based on 
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the assumption of a functional relationship between the 
obligors' expected default rates and different background 
factors. These background factors may represent either 
observable or unobservable variables. Observable vari- 
ables typically depend on the general state of the econ- 
omy, and the unobservable variables often represent some 
random risk factors [5D]. The event of default is often 
modeled by an intensity process, e.g. a Poisson process. 
Unlike the structural approach, the reduced form ap- 
proach is not directly based on a dynamical description 
of the economy. Reduced form models arc often imple- 
mented as black-box models, where the accuracy of the 
model outcome is more important than an intuitive eco- 
nomical interpretation of the mechanisms included in the 
black box. 

One of the key issues when setting up a credit risk 
model is to define an event that leads to an actual loss. 
Usually, credit risk modelers use cither of two definitions 
of credit loss: the default mode paradigm, or the mark- 
to-market paradigm. Within the default mode (DM) 
paradigm, a credit loss occurs only when an obligor de- 
faults within the bonds maturity time. Within the mark- 
to-market (MTM) paradigm a credit loss can occur with- 
out an actual default. The creditor can lose money when- 
ever an obligor's credit worthiness deteriorates. Both of 
these two modeling approaches are common in current 
vendor credit risk frameworks. Typical arguments for 
using one of them are the simplicity of the DM model 
and the multi-state nature of the MTM model. The DM 
model suits better for creditors which only are interested 
in a buy-and-hold portfolio, while the MTM model is 
more adequate for pricing decisions of more liquid cred- 
its [U]. 

There is a wide range of different credit risk manage- 
ment tools available in the financial industry. Many of 
those models seek to estimate the full distribution to be 
able to calculate different statistics measures [2D] . Exam- 
ples of vendor credit risk models are CreditMetrics (by 
RiskMetrics Group), PortfolioManager (by KMV), Cred- 
itPortfolioView (by McKinsey & Co), and CreditRisk + 
(by Credit Suisse Financial Products), see Ref. [22] and 
a review in Ref. 1231 . 



III. A MODEL WITH JUMPS AND 
CORRELATIONS 



We set up our structural credit risk model in Sec. Ill A 



In Sec. |IIIB[ we discuss its generality by relating it to 
a few other scenarios where structural models can find 
application. 



A. Setup of the model 

We model the time evolution of the asset value of every 
company by a stochastic differential equation of the form 



dV 



= /idt + aeVdt + dJ . 



(1) 



Apart from the jump term dJ, this is a geometric Brow- 
nian motion with a deterministic term [idt describing 
the exponential growth of the asset value and a stochas- 
tic term dW — eVdi representing the fluctuations as a 
Wiener process. Here, /j, is the drift, a the volatility (con- 
stant) and e an independently distributed random num- 
ber in each time step. We add the jump term dJ which 
is not contained in Merton's original model. In previ- 
ous works, jump-diffusion models have been considered 
for stock returns (see e.g. Ref. [23]) and also for credit 
risk [13, 25 26 , [27] . The economical interpretation of the 
jump term is that a major setback of an asset value is 
possible at any time. These setbacks are larger than the 
volatility admits and may be explained by events labeled 
as crises, originating from legal, operational, political or 
other external or internal factors. The probability that 
an economical setback takes place during the lifetime of 
a bond is typically very small. Well known examples of 
such economical setbacks for a large group of companies 
are the great stock market crash of 1929, the oil crisis 
of the mid-seventies and the "Black Monday" crash in 
1987. We model the jumps by a Poisson process with 
intensity A. We recall that, in such a process, the prob- 
ability function for the event to occur n times between 
zero and the time t is given by 



Poissoii 



Pn 



(*) 



(At)' 



exp (—At) 



(2) 



The size A of the jump, measured in units of the cur- 
rent asset value V(t), is a random variable with a distri- 
bution which we have to specify. Jumps can be positive 
or negative. The largest possible negative jump is 100% 
of the current asset value. Based on this information, 
a possible distribution of the jump size A is a shifted 
lognormal distribution, A + 1 ~ LN(/ij + 1,<7j), with 
mean fj,j and standard deviation <jj. A time series for 
the asset value including a negative jump is shown in 
Fig. [3] Without the jump term, the distribution of the 
asset price V(t) is log-normal. The jumps render the 
tails of the asset price distribution fatter. Fat tails are 
empirically observed [Tj. As this clearly affects the loss 
distribution, we find it important to include such jumps. 
The parameters of the jump process can be adjusted in 
order to match the tail behavior of a given empirical time 
series of the asset value. 

To determine if default has occurred, we compare the 
asset value at maturity T to the company's financial 
obligations, i.e. the face value F of the bond. This is 
in distinction to the so-called first-passage models, see 
e.g. Refs. pT5l 1231 12"8"] . where default occurs as soon as 
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FIG. 3: Two time series for the asset value V(t) versus time 
t. They are equivalent except that one of them includes a 
large negative jump given by the jump term J(t). The lower 
time series is damped because the factor V in Eq. |T| is much 
smaller after the jump. 



the asset value V(t), with < t < T, falls below some 
specified threshold function D (t). For simplicity, we do 
not choose such a function and work with the face value 
F. In Fig. |4j the asset value process with the thresh- 
old is visualized. We employ the default mode paradigm, 
which means that an actual loss only occurs in case of a 
default. The size of the loss in case of default is given by 
the face value of the loan minus the total asset value at 
maturity. We normalize the loss to the face value F, 



L 



F - V{T) 
F 



(3) 



implying that we have < L < 1 for the normalized loss 
L. Many studies assume, for the sake of simplicity, that 
everything is lost in case of default, i.e. L can only take 
the values and 1, see e.g. Refs. [28j|29]. In this simple 
case, many distinctive features of the loss distribution, 
as seen e.g. in Fig. [12J cannot be reproduced, and, par- 
ticularly, the tail behavior of the loss distribution is dif- 
ferent. Zhou [13] studied a first-passage model for credit 
risk based on a jump-diffusion process. This leads to a 
recovery rate which depends solely on the jump process, 
whereas in our model it is V(T)/F, i.e. it is determined 
by the whole asset process. In Ref. [13] default prob- 
ability and credit spread are studied in dependence of 
maturity T. The present study is different, because we 
we extensively examine the full loss distribution of the 
credit portfolio. 

It is important to include correlations between the 
companies in our model. We mention that financial cor- 
relations presently find considerable interest in the econo- 
physics community. These correlations are noise dressed 
if measured for finite time series. The impact of noise 
dressing on portfolio optimization and methods how to 
reduce the noise to find the true correlations are much 
discussed in the literature [30][3D[32][MllSl[35l[36l[3ll 
155] , The influence of noise dressing on credit risk has 
very recently been studied in Ref. [7] . To model the cor- 



relations, we employ Noh's model [39] which belongs to 
the class of factor and arbitrage pricing models [2H1 HQ] . 
Noh's model produces correlated normalized time series 
Mfe(t). We use the conventions of Ref. [35]. 

To achieve realistic portfolio correlations we divide the 
companies into different branches b = 1,2,. ..,£>. The 
total number of branches represented in a portfolio is 
given by the number B. The different companies in the 
portfolio are given by the index k = 1,2, ... ,K and the 
branch index & is a function of the company index, i.e. b = 
b(k). The number of companies in a specific branch b is 
given by Kb- For the companies that are in no branch we 
have 6 = 0, and the number of those companies is given 
by k. This means that we can write the total number of 
obligors K in a portfolio as 



K 



b=l 



(4) 



Within a branch b the Kb companies are assumed to be 
correlated with a specified correlation coefficient To 
achieve this, the one-factor model adds a part of the 
branch specific time series r]b(k){t) to the branch inde- 
pendent time series e k (t) . A sum of these time series is 
used to construct the asset returns for the companies in 
the portfolio. The correlated time series Mfc(i) has the 
form 



M k (t) = 



Pb(k) 



1 



- Pb(k) 




e k (t). (5) 



The entries in both ijb(k)(t) and are uncorrelated 

and standard normal distributed. The weights Pb(k) mea- 
sure the correlation and satisfy Pb(k) > 0. In particular, 
we have po = for those companies which are in no 
branch. The normalized, correlated time series Mfc(i) 
can be arranged in a K x T matrix M, where K is the 
number of obligors and T the length of the time series. 
The corresponding K x K correlation matrix C(T) is de- 
fined as 

C(T) = ijlfAft = (M k (t)Mi(t)) T , (6) 

where is the transpose of M. The index T on the 
brackets indicates that the average depends on the length 
T of the time series. If the time series are infinitely long, 
i.e. T — > oo, the correlation coefficient Cjy(oo) for com- 
pany k and I, 



C ki ((X)) = 



1 



1 



-Pb(k) 



(Pb(k)$b(k)b(i) + hi) ■ (7) 



This value should resemble the best estimate for the cor- 
relation and is referred to as the true correlation. The 
correlation matrix C(oo) consists of B square blocks of 
dimension Kb x Kb at the diagonal. The off-diagonal el- 
ements are Cj, = Pb/(1 + Pb) for branch b, and the no- 
branch companies are represented by an identity matrix. 
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TABLE I: Input/output of the model. 



Input 


Description 


Unit 


K 


Number of obligors in the portfolio. 


— 


T 


Timespan of the bond (maturity). 


Year 




Asset drift. 


[Year]" 1 


(7 


Asset volatility. 


[Year]" 172 


A 


Jump intensity. 


[Year] -1 




Mean of the jump size. 


Percent 


cj 


Standard deviation of the jump size. 


Percent 


F 


Face value of the bond. 


$ 


V 


Asset value at the issue date t = 0. 


$ 


C 


Correlation matrix. 




Output 


Description 


Unit 


Pd 


Probability of default. 


Percent 


P(L) 


loss distribution L. 





All diagonal elements of C(oo) are equal to one, and all 
entries that are not mentioned above are zero. However, 
for finite length T < oo, the true correlations are noise- 
dressed, because every matrix element carries a random 
number as offset and the block structure is obscured. 

Implementing Noh's model, we are able to analyze a 
wide range of different portfolio compositions. We em- 
phasize that we automatically include noise, because we 
work with time series of finite lengths. For this investiga- 
tion, we find this more realistic because it seems to match 
the present usage by practitioners. We notice that linear 
correlations between the asset processes are the simplest 
form to model the correlations between defaults. They do 
not account for credit risk contagion across firms and pe- 
riods of default clustering, see e.g. Refs. [51 151 IT21 I5T1 W2\ . 

Tab. [I] summarizes the model parameters and Fig. [4] 
shows a visualization of the underlying asset value pro- 
cess. The model comprises a large number of parameters. 
However, this is indispensable if one aims at setting up a 
realistic model. Credit risk is a problem that involves a 
high degree of complexity. Importantly, none of our pa- 
rameters is a hidden parameter. All of them have a direct 
interpretation and are measurable observables. Thus it 
is certainly possible to calibrate our model by fitting all 
parameters to real market data. Admittedly, this task 
might be time-consuming or even difficult, but it is def- 
initely feasible. A key purpose of the present study is 
to investigate the tail behavior of the loss distribution. 
Thus, even if it is difficult to determine some of the pa- 
rameters sharply, that is, only within some uncertainties, 
our model yields detailed information on how the tail de- 
pends on these parameters. 

B. Related scenarios 

The task of finding the probability to hit a certain 
threshold for an object, a particle, say, whose motion is 



V(t) 



V(0) 
F 




T t 



FIG. 4: A visualization of the underlying asset value process, 
involving the time series of the asset value V(t) versus time 
t with initial value V(0) = Vo, face value F and maturity 
T. The curve to the right is the distribution of the asset 
value at maturity. The shaded area corresponds to the default 
probability. 

described by a stochastic process closely relates to the 
much discussed problem of the stopping time distribu- 
tion. The stopping time is the time at which the parti- 
cle hits the threshold. The corresponding distributions 
are non-trivial objects, because they tend to have non- 
generic features [43] • This explains the considerable in- 
terest which they attract in mathematical statistics. 

An application in physico-chemistry is the following: 
A (Brownian) particle moves stochastically in a suspen- 
sion confined by a container. It is absorbed when it even- 
tually hits the wall of the container, or it reacts there, or 
it penetrates the wall if the latter is semipermeable and 
separates two suspensions. Obviously, the associated loss 
distribution is very similar to the one in credit risk, if 
we consider a first passage model. The numerical sim- 
ulations can be easily adjusted to such a scenario. One 
could even imagine a force acting between the stochastic 
particles which would imply a certain correlation between 
them. However, we do not elaborate on this further, be- 
cause the purpose of the present discussion is only to 
establish a qualitative connection. 

Another example is a farm where different types of 
crops are grown. The time evolution of their quality and 
their future value has a deterministic and a stochastic 
component. The analogue of default for a given type of 
crop arises when its value falls below a certain thresh- 
old. This could be due to bad weather, diseases, pests, 
fire, etc., also motivating the inclusion of jump processes. 
The time series for the crops are correlated, for example, 
because the diseases are contagious. The knowledge of a 
distribution for the loss similar to the one in credit risk 
would be most helpful for the farmer. 

The last example shows that the structural model for 
credit risk can find application in a large number of logis- 
tic problems where various objects or quantities have to 
be available in a certain state at some time in the future. 
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IV. ANALYTICAL DISCUSSION 



We investigate a simplified version of our credit risk 
model analytically. The purposes of this discussion are, 
first, to illustrate the general mechanisms and, second, 
to look at the tail behavior. In doing so we demonstrate 
how and with which speed the loss distribution converges, 
under certain assumptions, to a universal limit as the 
number K of companies is made large. Thus, in this case 
study, we do not include the jump term in the diffusion 
process. Moreover, we also make the assumption that 
the obligors are uncorrelated. We distinguish between 
individual and portfolio losses in Sees. |IV A| and |IVB| 



A. Individual losses 

The asset value V(t) for every company k with k — 
1, . . . , K follows a geometric Brownian motion; its distri- 
bution is log-normal. With the initial value V(0) = Vq at 
t = 0, the distribution of V(T) at maturity t = T reads 



PkiViT)) = 



1 



exp 



V(T)V2wa 2 T 

(HV(T)/V )-^-ay2)TY 
2a 2 T 



(8) 



We use the company index only for the distribution p^, 
and we suppress it in the quantities of the asset value 
process. Due to < V(T) < F at maturity in the case 
of default, we have to truncate Eq. pi. Using Eq. ([3| 
we map the distribution of V(T) to that of the loss given 
default L for this company, 



Pk(L) 



exp 



1 



P D (1 - L)V2tt<t 2 T 
(ln(F(l-L)/y )-(^-a 2 /2)T)^ 
2a 2 T 



The factor needed to restore the normalization 

11 fHF/V )-^~a 2 /2)T 
D 2 2 



V2^T 



.(9) 



(10) 



is the default probability for this company fc, where we 
employ the error function according to the definition 



erf (x) = -j= J exp(-£ 2 R 



(11) 



The number is the probability that the asset process 
V(t) is below the threshold at maturity T, i.e. for having 
V(T) < F. The default probability is shown in Fig. [4] as 
shaded area. 

The impact of the drift term (i and the volatility a on 
the default probability is shown in Fig. ([5j. The range 



FIG. 5: The default probability ( 10 1 as a function of the 
parameters fi and a for the values Vb = 100, F — 75 and 
T — 1. The default probability tends to increase rapidly 
when choosing parameter values that are too large. Also one 
can see that it is an even function in a, indicating that it is 
only the absolute value that matters. 



in which realistic default probabilities are generated is 
vary narrow, both in the /j, and a direction. For our 
simulations we choose \x = 0.05 and a = 0.15, which 
results in a default probability of Pd ~ 0.0148. 
For later purposes, we calculate the n-th moment 



(L n k ) = J L n Pk (L)dL 



(12) 



of the distribution ([9]). By using the index k in our nota- 
tion, we want to underline that the (L^) are the moments 
of the distribution for the individual loss given default. 
We arrive at 



U 3=0 J 



a 2 ( F 

exp ( j(j - l)yT + j {[iT - In — 



1 + erf 



(l-2j)a 2 T/2-^T + ln(F/V ) 



(13) 



V2W 2 T 

after a straightforward calculation. 

B. Portfolio losses 



The loss L of the total portfolio is the arithmetic mean 
of the individual losses, 



1 K 
L= xYl LkIk 



(14) 



k=l 



where is the loss given default for the individual bond 
with index k and K is the total number of bonds in the 
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portfolio. As we are interested in the distribution of the 
portfolio loss, we have to introduce the default indicator 
I k for company k by 



h = 



if 
if 



V(T) < F 
V(T) > F 



(default) 
(no default) 



(15) 



As the distribution of the default indicator for company 
k, we choose 



Pk{h) - (1 - PnWk) + P D S(h - 1) 



(16) 



The default probability ( |10| depends on the specific pa- 
rameters which go into the asset value process for com- 
pany k. To obtain the distribution of the portfolio loss 
we have to average over all distributions of the individ- 
ual losses p k (L k ), given by Eq. ([9]), and over the indicator 
distributions (|T6l, 



+ oo 



P(L) 



J dhpiih)--- J dI K p K {I K ) 

so — oo 

1 1 
J dLipi(Li) ■ ■ ■ J dL K p K {L K ) 



*( L -^X>^-) ■ ( 17 ) 



We notice the subtle difference between the loss distribu- 
tion and the distribution of the loss given default. This 
difference is best understood by considering Eq. (17 1 for 
K — 1 , which yields 



p{L) = (1 - P D )6(L) + P DPl (L) 



(18) 



This is the weighted sum of the distributions for the case 
of no default and for the case of default. 

We calculate an asymptotic approximation to p{L) for 
a large number K of companies, i.e. for large portfolio 
size. Here, we make the further assumptions that the 
face values and the parameters of the geometric Brow- 
nian motion are the same for all companies k. These 
additional assumptions are not strictly necessary for the 
feasibility of the calculation, but they make the resulting 
expressions very compact. Details of the computation 
are given in Appendix [A] We obtain 



P(L) 



1 

2n 



dcuexp (—icu (L — Pjj(L^))) 



'■■■■■V\-^(PD(L 2 k )-P 2 D (L k ) 



exp 



6K 2 



{Pd(L 



-3P 2 (L k )(L 2 ) + 2P 3 D (L k 



(19) 



as asymptotic approximation to order 1 / K 2 . If we skip 
the 1 / K 2 term and settle with an 1/K expansion, we find 
the shifted Gaussian 



P(L) 



K 



exp 



2n(P D (L 2 ) - P 2 (L k y) 
K(L-P D {L k )f 



2(P D (Ll) - P 2 (L k y) 



(20) 



Thus, the Central Limit Theorem applies if the number 
K of companies is very large. The expected loss EL, 
i.e. the peak position, is the mean value of L kl weighted 
with the default probability Pjj. This is so, because we 



defined the portfolio loss L as the arithmetic mean ( 14 1 



The unexpected loss UL quantifying the risk is given by 
the square root of (Pd(L 2 ) — Pp(L k } 2 ) / K which becomes 
smaller with K. Even in the Gaussian limit, the depen- 
dence of p{L) on the parameters F, Vq, /i, a and T is 



non-trivial due to the rather involved expressions ( 10 1 



for the default probability and ( 13 1 for the moments. We 
will return to this point. Due to the specific nature of 
the approximation that leads to Eq. ( 19 1 , the first mo- 



ments of the exact distribution are preserved up to the 
highest order of lo which is considered in the exponential. 
Therefore, Eq. (19 1 contains the correct first, second and 
third moment, while the Gaussian limit in Eq. (20 1 still 



has the correct mean and variance. 

Further, it is important to note that expression (19 1 



holds for any structural model for uncorrelated and ho- 
mogenous portfolios, i.e. with all face values and param- 
eters of the geometric Brownian motion being the same 
for all companies k.. It does not depend on the choice of 
random processes for the asset values V(t), and it is suf- 
ficient to know the default probability and the first three 
moments of the individual loss given default distribution. 
It is straightforward to generalize Eq. (19 1 to the case of 
inhomogenous portfolios (see Appendix |B[). 



C. Exact loss distribution versus approximations 

The approximation in the previous subsection has the 
advantage of being easily extendable to more general 
portfolios and exactly conserving the first moments of 
the distribution up to the order of approximation. How- 
ever, it does not reproduce the shape and, in particular, 
the tail behavior very well. 

A better approximation of the loss distribution is pos- 
sible, if we restrict ourselves again to homogenous port- 
folios. Using the default indicator function defined in 
Eq. (16 1, we can then rewrite Eq. (17 1 as a combinato- 



rial sum 



K 

E 

3=0 



(1 



Pn) K ~ j P 3 D Fj(L) 



(21) 
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where we define the function Fj(L) as 

1 1 
Fj(L) = J dLwil*)'-- J 'II.jPJLj 
o o 

« y~ I duexp (-iu (l - j^(L k ) 



(22) 



-iM-|^((4)-(^) 2 ) 



6K 3 



((Ll)-3(L k )(Ll) + 2(L k ) i 



The approximation for Fj(L) follows the same line of 
arguing which is outlined in Appendix [A| Note that each 
term of the sum in Eq. pTj ) corresponds to the event 
that exactly j defaults occur. In particular, for j = 
the delta peak (1 — Pd) K 8(L) is obtained exactly in this 
approximation . 

In Figure [6] we compare the exact loss distributions 
for K = 10,100 and 1000 to the asymptotic approxi- 
mation ( (T9| ) to order 1/K 2 , and to Eq. (211 with the 
approximation in Eq. (p2|. In both cases, the analytical 



approximations arc evaluated numerically The exact loss 
distributions are obtained by calculating numerically the 
convolutions in Eq. ( 22 ) , using the characteristic func- 
tion of Pk(L), and inserting these Fj(L) into Eq. (21 1. 
As numerical values for the model parameters we choose 
fi = 0.05, a = 0.15, T = 1, V = 100 and F = 75. 

For K = 10 the delta peak (1 - P D ) K 5{L) (not shown 
in the plots) is quite dominant and leads to a very poor 
result of approximation (19 1, while the approximation of 
Eq. (21 1 already yields a quite reasonable result. For K — 
100, features of the Gaussian limit are already present. 
Pictorially speaking, the Gaussian moves from the left 
into the picture. The approximation of Eq. ( |21| captures 
the tail behavior of the distribution quite nicely. 

For K = 1000, the Gaussian limit (20 1 is almost 
reached, but, interestingly, the distribution is still slightly 
asymmetric. The agreement with the asymptotic approx- 



imation (U9 



is convincing, although the tail behavior is 



not as well described as by the approximation of Eq. (21 1. 

Figure [6] yields a good intuition for the tail behavior 
and for the speed of convergence to the Gaussian, or any 
other universal limit. Remarkably, even for K = 1000 
slight deviations from the Gaussian are seen. This means 
that, importantly, only really large portfolio sizes imply 
universal shapes. 

A measure for the tail behavior of the loss distribution 
is the kurtosis, which is defined as a normalized form of 
the forth central moment p,<± 



ft - ^ 

\4 



(23) 
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FIG. 6: Loss distribution for three different portfolio sizes 
K = 10, 100, 1000, respectively. The insets show logarith- 
mic plots. Black solid lines show the analytical resu lts, while 
the blue dashed-dotted lines correspond to Eq. (191, and the 
orange dashed lines to the approximation of Eq. (21 1. The 
model parameters are fi = 0.05, a = 0.15, T = 1, Vo = 100 
and F = 75. 



where [12 is the second central moment. Since the kur- 
tosis of the normal distribution is equal to three, the 
kurtosis excess is defined as 



72 



= ^-3. 



(24) 



It can be shown analytically that the kurtosis excess for 
homogenous and uncorrelated portfolios is proportional 
to 1 / K. Such considerations are important in modern 
portfolio theory, for example, when one tries to minimize 
risk by diversification. This can require an enlargement 
of the portfolio. 
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We notice that the symmetric shape of the Gaussian 
can only be reached without correlations between the as- 
set processes: with correlations, the probability of default 
is enlarged without compensation on the positive side, 
because in the best case, none of the obligors defaults. 
Thus, the distribution is asymmetric. 



For example, the two first moments are, for n = 1 

A Fi R ) \\ J/K-l 



-<K) 



(30) 



D. Drill down risk 

Consider a portfolio consisting of K companies. The 
moments of the loss given default distribution are as- 
sumed to be known. In which way will the properties of 
the portfolio change if one decides to add or remove a 
company? If one had this information an optimal strat- 
egy could be developed to drill down the parts within the 
portfolio which are the most risky, or to decide which new 
company should be added, so that the additional risk is 
minimal. 

Since we want to allow the individual face values Fk 
and the initial values Vok to be different now, it would 



be unsuitable to use the dimensionless definition (14) of 
the losses. Instead we use a more general definition of 
the individual loss 



Tfc = LkFk — Fk 

which now has the dimension of dollars 
for the entire portfolio of size K is then 

'A' 



Vfe(T), (25) 
The total loss 



(26) 



which again is a normalized quantity with < L^ K ' < 1. 
Ik is the indicator function for company k, as defined in 
Eq. ( 15 1. If all face values Fk are the same, Fk = F, the 



loss will again reduce to Eq. (|14j, since 



Hk=l T khk _ FkY,k=l L k I k _ 1 



A' 



L (K) 



EfcLi F k 



KF k 



A 
k=l 



LkFk 



(27) 

Now, is it possible to express the moments of the loss 
distribution for K companies, in terms of the correspond- 
ing loss distribution for K — 1 companies? Indeed it is. 
See Appendix [C] for details. The n:th moment can be 
expressed as: 



Ml 



A 




A-l 



(Ik) TO 
(p{K-i)y 



(28) 



where v is an integer and F^ K ^ is the sum over all the K 
face values: 



F {K) 



K 



(29) 



and for n = 2 



A 



F {K) 



A-l 



2 ^(/A)(r,)((^)>A- 1 



(F( 



1 



(F( 



{Ik){T%). 



(31) 



V. NUMERICAL DISCUSSION 

The numerical analysis of our model is done with 
Monte Carlo simulations. To achieve results that are 
statistically reliable, we simulated between 10 4 and 10 5 
asset processes. An important remark is in order. Given 
the rich variety of credit contracts, we decided not to 
calibrate our model with real market data, because this 
would match only one particular scenario. However, we 
have chosen the parameter values in such a way |16j that 
the model outcome is economically realistic. In particu- 
lar, the parameters of the jump process have been cho- 
sen in order to reproduce a realistic tail behavior. For 
homogenous and uncorrelated portfolios the loss distri- 
butions obtained by the Monte Carlo simulations agree 
very well with the analytical expectation. In Sec. |V A| we 
study the impact of maturity, drift and volatility. Lever- 
age, jumps and correlations are addressed in Sees 



VB 



|V C| and |VD[ respectively. We discuss the competition 
between jumps and correlations in Sec. |VE] 



A. Maturity, drift and volatility 

To give the reader a better feeling for the sensitivity 
of the loss distribution, we first take a look at the model 
in its simplest form, i.e without jumps and correlations, 
and study its dependence on maturity T, drift /i and 
volatility a. This is again done by calculating Eq. (21 1 



numerically, while in the remainder of the paper we use 
Monte Carlo simulations. 

We begin with the maturity T and keep fi = 0.05, a = 
0.15, Vo = 100 and F = 75. It is quite common that loans 
span over a period longer than T — 1 year which has been 
considered so far. In Fig. [7] the expected loss EL and the 
unexpected loss UL are shown as a function of maturity 
T. The expected loss EL is the mean value Po{Lk) and 
the unexpected loss UL squared is the variance (Pd (-kfe) - 
Pp(Lk) 2 )/K, where the default probability Pd and the 
moments (LJ?) are given by Eqs. (10 1 and (13). Closer 
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FIG. 8: Expected loss EL and unexpected loss UL as a func- 
tion of the drift fi and the volatility a for K — 1000. 



FIG. 7: Expected loss EL and unexpected loss UL as a func- 
tion of the maturity T for K = 1000. 



inspection shows that EL grows monotonically in T to its 
saturation value unity, if, for fixed (positive) drift \x and 
volatility a, the condition fi < a 2 /2 is met. However, if 
/i > er 2 /2, EL has a maximum for some finite value of T 
and the saturation value of EL is zero for T — > oo. For the 
parameter choices fj, = 0.05 and a — 0.15, the maximum 
is at T « 12.56, as seen in Fig. [7] Similar considerations 
apply to UL, with the maximum at T w 17.55. 

In the sequel, we put T = 1 and investigate the de- 
pendence on drift /i and volatility a. We choose 0.05 < 
H < 0.15 and 0.15 < a < 0.35 which is motivated by 
economical data. The expected loss EL and the unex- 
pected loss UL for K = 1000 are depicted in Fig. [8] For 
K = 1 (not shown), both EL and UL are proportional 
to the volatility and only weakly dependent on the drift. 
For portfolio sizes K > 1 and fixed volatility a, both EL 
and UL increase with decreasing \i. This reflects that 
the probability to default increases with smaller \x. Im- 
portantly, the loss distribution becomes more sensitive to 
this effect with growing portfolio size. This is so, because 
the portfolio loss is an average of individual losses. The 
expected loss EL and the unexpected loss UL increase 
with volatility. As for the drift, the sensitivity of the 
loss distribution to the volatility grows with the portfo- 
lio size. A larger volatility implies larger asset returns 
which cause a greater activity of the asset value. This, 
in turn, produces excessive losses. 

We conclude that the qualitative behavior of the loss 
distribution as function of drift and volatility is very sim- 
ilar. Economically, this means that a downwards change 




0.05 0.10 
drift 



FIG. 9: The kurtosis excess as a function of the asset drift 
(j, for portfolio sizes K — 10 (dotted), K — 100 (dashed) and 
K = 1000 (dashed-dotted). 



in the general trend of the asset prices during a quiet, 
not volatile period affects the credit markets in a sim- 
ilar way as an upwards change in the market activity, 
i.e. volatility, in times of a stable general trend in the 
markets. According to Eq. (10 1, the default probability 



is almost proportional to both drift and volatility in the 
parameter range chosen. 

Finally, we study how drift \i and volatility a affect 
the tail behavior of the loss distribution. We recall that 
the universal limit is in the present case always Gaussian 
for very large K. It can be shown analytically that the 
kurtosis excess of uncorrelated portfolios scales as 1 /K. 
Figures [9] and [10] show the kurtosis excess as functions of 
drift and volatility, respectively. In Fig. [9] one sees that 
a higher drift gives loss distributions with fatter tails. 
This is more pronounced for large portfolio sizes. Fig- 
ure [10] shows that the kurtosis excess approaches zero, 
i.e. the distribution becomes more Gaussian, with grow- 
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FIG. 10: The kurtosis excess as a function of the asset volatil- 
ity a for portfolio sizes K = 10 (dotted), K — 100 (dashed) 
and K = 1000 (dashed-dotted). 



ing volatility. Moreover, one also observes that the speed 
of convergence to a Gaussian is higher for larger portfolio 
size, which reflects the Central Limit Theorem. 



B. Leverage 

The ratio F/Vq between the face value of the bond F 
and the initial asset value Vq is referred to as leverage. 
It is very important for credit risk managers to know 
its influence on the loss distribution. In our model we 
are able to set both the face value of the bond and the 
initial asset value. However, as it suffices to study only 
the leverage, we keep Vb = 100 fixed and only vary the 
face value F. The portfolio default probability PD, ex- 
pected loss EL and the unexpected loss UL grow with 
the leverage. This is so, because the higher the leverage, 
the more likely are early defaults. It is a remarkable re- 
sult that large portfolios are extremely sensitive to the 
leverage. As every realistic portfolio contains bonds with 
different leverage, it is of considerable interest to study 
how the loss distribution changes if the leverages are dis- 
tributed. Here, the (relative) loss L is defined as the 
sum of the absolute losses FkLklk divided by the sum of 
the face values Fk of the companies k = 1, . . . , K. Here 
Lk is the (relative) loss given default for company k and 
Ik the corresponding default indicator. For convenience, 
we choose uniform distributions with width AF centered 
around F = 75. The results are shown in Fig. 11 Ex- 



pected loss EL and unexpected loss UL grow with AF. 



C. Jumps 

The need to include jumps in a realistic model has 
been pointed out in Sec. |III A| Jumps are important de- 
spite the fact that the jump intensity is typically very 
small. Reasonable economical values of the jump inten- 
sity are around 0.01 jumps per year. If not otherwise 
stated, we use the following values for the jump process: 
jump intensity A = 0.01, mean jump size /ij = —0.4, 
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FIG. 11: Loss distributions for portfolios in which the lever- 
age is uniformly distributed. The face values are random in 
a window of width AF around F — 75 with AF — (solid 
lines), AF = 10 (dashed lines) and AF = 20 (dotted lines), 
for three different portfolio sizes K — 10, 100, 1000, respec- 
tively. The expected loss EL is 0.076%, 0.095% and 0.157% 
for the three different leverages (independent of K) . The in- 
sets show, for every portfolio size, the portfolio default prob- 
ability PD and the unexpected loss UL from top to bottom 
for AF = 0, 10, 20. 



and standard deviation oj — 0.3. In Fig. [12] we plot 
the loss distributions for three different jump intensities, 
A = 0.005,0.01,0.015. Not surprisingly, PD, EL and UL 
increase with the jump intensity. As for any uncorre- 
lated portfolio, EL is independent of the portfolio size, 
while UL scales with 1/y/K. For K = 1000, PD has 



already saturated at 100%. In Figure 13 the kurtosis ex- 



cess is depicted as a function of the jump intensity A. 
The kurtosis excess grows with the jump intensity up to 
a certain value of the intensity, then it decreases again 
and approaches zero asymptotically. For larger values of 
A, the loss distribution is more smeared out, because the 
losses produced by the jumps dominate the defaults due 
to the standard geometric Brownian motion. In other 
words, the emergence of the maximum can be explained 
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FIG. 12: Loss distributions for portfolios with different jump 
intensities A = 0.005 (solid lines), A = 0.01 (dashed lines) 
and A = 0.015 (dotted lines) for four different portfolio sizes 
K — 10, 100, 1000, respectively. The expected loss EL is 
0.15%, 0.22% and 0.29% for the three different jump inten- 
sities. The insets show, for every portfolio size, the portfolio 
default probability PD and the unexpected loss UL from top 
to bottom for A = 0.005, 0.01, 0.015. 
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FIG. 13: The kurtosis excess as a function of the jump in- 
tensity A. The result is shown for K = 10, and it scales with 
l/K. 



as follows. The kurtosis is high for a distribution that has 
a sharp peak around the mean and a fat tail. When the 
intensity A grows starting from A = 0, the tail becomes 
fatter, because the kurtosis increases. But at some value 
of A, the tail becomes so heavy that it interferes with 
and ultimately destroys the sharp peak of the distribu- 
tion. Hence, the kurtosis excess must have a maximum 
value for some intensity A. Furthermore, we notice that 
the kurtosis excess scales as l/K with the portfolio size. 

As discussed in Sec. |III A] we make the specific assump- 
tion that the distribution of the jump sizes is log-normal. 
Hence, it is important to see how the loss distributions 
depend on the parameters of the jump size distribution. 
The losses grow when the absolute value of the jump 



mean fj,j increases. This is illustrated in Fig. 14 As for 
the jump intensity, the sensitivity of the loss distribution 
to the jump mean does not significantly depend on the 
portfolio size. Furthermore, as displayed in Fig. 15 the 



credit losses slightly increase, for all portfolio sizes, with 
the jump size standard deviation oj. 



D. Correlations 

To begin with, we choose the simplest possible correla- 
tion matrix C = C(T), i.e. a portfolio where all com- 
panies are in the same branch. Figure [16] shows the 
structure of the correlation matrix and the loss distribu- 
tions for different portfolio sizes and branch correlations 
c = C'b- In the sequel, the structure of the correlation 
matrix is always presented as a gray scale image where 
the intensity corresponds to the correlation. Black areas 
stand for unit correlation, white areas for zero correla- 
tion. The expected loss EL is 0.076% for all three branch 
correlations, i.e. it is independent of the branch correla- 
tion and, as in the case of uncorrelated portfolios, it is 
also independent of the portfolio size. The unexpected 
loss UL becomes larger as the branch correlation grows, 
while the portfolio default probability PD decreases. The 
latter can be understood, because the portfolio increas- 
ingly acts as a single company with growing correlation 
strength. The portfolio default probability makes a tran- 
sition from 1 — (1 — Pd) K , corresponding to K uncor- 
related companies (c — 0), to Pjj for a single company 
(c=l). 

Figure [L6| shows that, as for a majority of the model 
parameters, the impact of the correlation on the loss 
distribution increases with the portfolio size. Remark- 
ably, there is a sizable deviation from a Gaussian-type-of 
shape. This is seen in Fig. [17] which shows the kurtosis 
excess of the loss distribution as a function of the branch 
correlation c — Cb for this one-branch scenario. The 
kurtosis excess for a portfolio of size K shows a transi- 
tion from the value 264.6/if for an uncorrelated portfolio 
(c = 0) to 264.6, which corresponds to the value for an 
individual obligor. Between c = 0.6 and c — 1 the kurto- 
sis excess even exceeds this value. To gain further insight 
into the influence of correlations, we now assume that the 
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FIG. 14: Loss distributions for portfolios with different mean 
values of the jump size, fj,j = —0.3 (solid lines), fj,j — —0.4 
(dashed lines) and \ij = —0.5 (dotted lines) for four different 
portfolio sizes K = 10, 100, 1000, respectively. The standard 
deviations is always erj = 0.3. The expected loss EL is 0.18%, 
0.22% and 0.26%, and the portion of negative jumps are 86%, 
91% and 94% for the three fij values, respectively. The insets 
show, for every portfolio size, the portfolio default probability 
PD and the unexpected loss UL from top to bottom for fij = 
-0.3,-0.4,-0.5. 



FIG. 15: Loss distributions for portfolios with different stan- 
dard deviations oj of the jump size, oj = 0.2 (solid lines), 
a j = 0.3 (dashed lines) and aj = 0.4 (dotted lines) for three 
different portfolio sizes K — 10, 100, 1000, respectively. The 
mean is always fij — —0.4. The expected loss EL is 0.20%, 
0.22% and 0.24%, and the portion of negative jumps are 96%, 
91% and 87% for these aj. The insets show, for every portfolio 
size, the portfolio default probability PD and the unexpected 
loss UL from top to bottom for oj — 0.2, 0.3, 0.4. 



size Ki of the branch is smaller than the total size K of the 
correlation matrix. Thus, K\ companies are correlated, 
K — Ki are not. In Fig. 18 we plot loss distributions 
for different branch sizes «i, the branch correlation is 
c = Cb = 0.5. We observe that bigger branch size yields 
a higher unexpected loss UL and a lower portfolio default 
probability PD. Importantly, the loss distributions for 
K = 1000 show a transition from a strongly leptokurtic 
to a more Gaussian-type-of, but still asymmetric, distri- 
bution as Ki decreases. This clearly visualizes the effect 



already mentioned in Sec. IV A A symmetric Gaussian- 
type-of shape can only be reached without correlations 
between the asset processes. 

In a realistic economical scenario, companies will be- 
long to B different branches. In Fig. 19 results are de- 



picted for portfolios consisting of companies in B — 1, 2, 5 
branches, all with the same branch correlation c — Cb — 
0.5. The structure of the correlation matrices is also 
shown in the figure. The number of branches does not 
significantly affect the loss distribution. Moreover, the 
loss distribution does not approach the normal distribu- 
tion as the size of the portfolio increases. This is so, 
because the high branch correlation c = Cb — 0.5 makes 
the branches in the portfolio behave like individual com- 
panies. Hence, the K xK correlation matrix is effectively 
only B x B. The curves in the subfigures of Fig. fl9] are 
almost identical, since for high c = Cb there is no signifi- 
cant difference between portfolios containing one, two or 
five companies. 

An interesting question is whether a small, highly cor- 
related branch has more or less effect than a large branch 
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FIG. 16: Loss distributions for portfolios with different 
branch correlations. The top left plot shows the structure 
of the correlation matrix, all companies are in one branch, 
indicated in grey. The other three plots show the loss distri- 
butions for three different portfolio sizes K = 10, 100, 1000, 
respectively. The correlations are measured by the branch 
correlation parameter c = 0.2 (solid lines), c = 0.5 (dashed 
lines) and c = 0.8 (dotted lines). The expected loss is 0.076% 
for all three branch correlations, i.e. it is independent of the 
branch correlation. The insets show, for every portfolio size, 
the portfolio default probability PD and the unexpected loss 
UL from top to bottom for c = 0.2, 0.5, 0.8. 



with weakly correlated companies. To look into this, we 
investigate how the loss distribution depends on the prod- 
uct K\c of the branch size «i and the branch correlation 
c = Cb- Figure [20] shows loss distributions for portfolios 
where 90%, 50% and 10% of the obligors in the portfo- 
lio are in a branch. The correlations for these branches 
are 0.1, 0.18 and 0.9 respectively. Thus, the product is 
Kic = 9 in all three cases. We notice that K\ denotes the 
branch size in percent of the total number K of obligors. 
We conclude that a small, highly correlated branch and a 
larger branch with less correlated companies are equiva- 
lent from a creditor's point of view. The unexpected loss 
UL is very similar, although not identical, for the three 
different parameter settings, while the portfolio default 
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FIG. 17: The kurtosis excess as a function of the correlation 
parameter c for portfolio sizes K = 10 (solid line) and K = 
100 (dashed line). The horizontal dotted lines indicate the 
kurtosis excess of the corresponding uncorrelated portfolios, 
264. 6/K, for K = 1, 10 and 100, respectively. 



probability PD is decreasing with increasing n\. The ef- 
fect on PD is observed to be strongest for intermediate 
size portfolios (K = 100 in our case). 



E. Jumps versus correlations 

Here, we study the interplay between jumps and corre- 
lations. In Figs . [21] and [22] the portfolio default probabil- 
ity PD and the unexpected loss UL are displayed as func- 
tions of the jump intensity A and the branch correlation 
c. To keep the discussion transparent, the correlation 
structure is a single branch containing all companies. 
The portfolio default probability PD increases with the 
jump intensity and decreases with growing correlation. 
For fixed small jump intensity A both PD and UL show a 
transition from K uncorrelated obligors to an individual 
one. However, for larger values of A the influence of the 
uncorrelated jumps becomes more dominant leading to a 
less pronounced dependence on the correlation strength 
c. The unexpected loss UL increases both with the jump 
intensity A and with the correlation c. As the correlation 
changes from to 1, the dependence of UL on the jump 
intensity changes from a sqare-root-type-of behavior to 
a linear one. The expected loss EL (not shown) grows 
linearly with the jump intensity and is independent of 
the correlation strength and the portfolio size. 



F. Correlated jumps 

In the previous section, we considered a correlated dif- 
fusion process with uncorrelated jumps. Now we study 
the case where also the jumps are correlated, i.e. a jump 
can occur both in the branch specific time series f7f>(/c)(£) 
and in the branch independent time series The 
jump sizes are scaled with the correlation coefficient ac- 
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FIG. 18: Loss distributions for portfolios with different 
branch size. The top left plot shows the structure of the 
correlation matrix, only the companies within the grey block 
are in the branch. The other three plots show the loss distri- 
butions for three different portfolio sizes K = 10, 100, 1000, 
respectively. The branch sizes are (in percent of K) ki = 30% 
(solid lines), Ki = 60% (dashed lines) and ki = 100% (dot- 
ted lines). The expected loss is 0.076% for all three branch 
sizes. The insets show, for every portfolio size, the portfolio 
default probability PD and the unexpected loss UL from top 
to bottom for m = 30%, 60%, 100%. 



cording to Eq. (|5j). We consider a correlation structure 
with five branches of size K/5. In Figs. 23 and 24 the 



portfolio default probability PD and the unexpected loss 
UL are displayed as functions of the jump intensity A and 
the branch correlation c. As the correlation changes 
from to 1, both PD and UL show a transition from K 
uncorrelated obligors to the behavior of the five branches. 
For higher jump intensities, this transition can be non- 
monotonous — for K — 10 we observe a maximum in 
PD at c Pi 0.3 and a minimum in UL at c w 0.35. It is 
important to note that the unexpected loss UL increases 
tremendously with the correlation coefficient. The ex- 
pected loss EL (not shown) grows linearly with the jump 
intensity and is independent of the portfolio size. For 
fixed jump intensity it shows a minimum at c = 0.5 due 



FIG. 19: Loss distributions for portfolios with different num- 
ber B of branches. The left column shows the structure of 
the correlation matrices, labeled CI, C2, C3, respectively. 
The right column shows the loss distributions for three differ- 
ent portfolio sizes K = 10,100,1000 for CI (solid lines), C2 
(dashed lines) and C3 (dotted lines). The branch correlation 
is c = Cb = 0.5. The default probability is 1.49%. The insets 
show, for every portfolio size, the expected loss EL and the 
unexpected loss UL from top to bottom for CI, C2, C3. 



to the scaling of jump sizes according to Eq. ([5|. 

In Figure [25] we show how the loss distributions are 
affected by introducing correlations. We compare three 
cases: uncorrelated jump diffusion, correlated diffusion 
with uncorrelated jumps and fully correlated jump dif- 
fusion, as discussed above. As correlation structure we 
choose again five branches of size K/5, and the corre- 
lation coefficient for each branch is set to c = 0.5. In 
the case of correlated jumps, we rescaled the jump inten- 
sity and jump sizes, so that they match the uncorrelated 
case. The correlations in the diffusion already lead to 
fatter tails; introducing correlations between the jumps 
enhances this effect. For the case of only five branches, 
it is most visible in smaller and medium sized portfolios. 
Figure [26] shows the loss distribution for K = 1000 and 
a correlation structure with 50 equally sized branches. 
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FIG. 20: Loss distributions for portfolios where the product 
of the branch size and the branch correlation is constant. The 
top left plot shows the structure of the correlation matrix. 
The other three plots show the loss distributions for three 
different portfolio sizes K — 10, 100, 1000, respectively. The 
branch sizes are (in percent) ki = 10% (solid lines), ki = 50% 
(dashed lines) and K\ = 90% (dotted lines). The expected 
loss is 0.076% for all three branch sizes. The insets show, 
for every portfolio size, the portfolio default probability PD 
and the unexpected loss UL from top to bottom for n\ = 
10%, 50%, 90%. 



The correlation of the diffusion terms leads to an only 
slightly fatter tail of the distribution, while the correla- 
tion of the jumps has a much more pronounced effect on 
the tail behavior. 



G. Drill down risk 

Finally, we consider how different indicators behave 
when removing one company from a portfolio of initial 
size K. We start with a portfolio which is composed in a 
rather realistic way. It consists of a large portion of low 
risk bonds, but also contains a portion of more risky ones. 
The portfolio was divided into five different categories 




K=100 




K=1000 




FIG. 21: Portfolio default probability PD as a function of the 
jump intensity A and the correlation parameter c for different 
portfolio sizes K = 10, 100, 1000, respectively. 



C = 1,2,..., 5. Each category has its own setup of the 
initial values Vq and face values F. All companies within 
a single branch are identical. The portfolio composition is 
summarized in Tab. (JTT|) . The fraction of the K companies 
that belongs to a certain catagory is a and the fraction 
of the total amount of money invested in that category 
is denoted 7. 

Next it has to be decided which company should be 
removed from the portfolio. A ranking system is estab- 
lished with respect to the product of the default proba- 



18 




bility and the mean loss 

Rk = PokfJ'k 



(32) 



for company k. The company with the highest value of 
Rk is considered to be the worst, and is the one which 
is removed. By using this ranking system one of the 
companies in category five will be removed, since they 
have the highest default probability and expected loss in 
the entire portfolio. 

The different indicators can be seen in Fig. (27), as 
a function of K. Naturally, smaller portfolios are more 



sensitive to the pruning procedure. For example, when 
going from K — 50 to K = 49, the EL is reduced by 
about 16% and the UL is reduced by an enormous 18%! 
The same comparison when going from K — 1000 to K = 
999 is less dramatic but none the less quite significant. 
In this case the EL is reduced by 0.75% and the UL by 
0.88%. 

The shape of the loss distribution is also affected 
slightly when removing only one single company. Both 
the skewness and kurtosis are increased, when going from 
K to K — 1 companies. On the other hand, both the 
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FIG. 25: Loss distributions for uncorrelated jump diffusion 
(solid lines), correlated diffusion with uncorrelated jumps 
(dashed lines) and correlated jump diffusion (dotted lines) 
for three different portfolio sizes K — 10, 100, 1000, respec- 
tively. The insets show, for every portfolio size, the portfolio 
default probability PD and the unexpected loss UL from top 
to bottom for the uncorrelated case, correlated diffusion with 
jumps and correlated jump diffusion. 



FIG. 24: Unexpected loss UL as a function of the jump inten- 
sity A and the correlation parameter c for different portfolio 
sizes K = 1, 10, 100, 1000, respectively. 



99.9:th percentile cvo.999 and the economical capital EC 
are decreased. 



1. Drill down risk in the presence of correlations 

How are the moments, and other indicators changed 
when there are correlations involved? In this section the 
portfolio, composed by the rule set in Tab. [TTJ is divided 
into ten different branches. Each branch is correlated 
with correlation strength one half (Cb = 0.5). Note that 
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Vo 
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Q 


F/V 


7 


1 


75 


50 


0.5 


0.67 


0.355 


2 


100 


75 


0.3 


0.75 


0.319 
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125 


100 


0.1 


0.80 


0.142 
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150 


125 


0.08 


0.83 


0.142 
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175 


150 


0.02 


0.86 


0.042 



TABLE II: Portfolio distribution, where F is the face value 
and Vo is the initial value of the company. The fraction of 
the K companies that belongs to a certain catagory is a and 
the fraction of the total amount of money invested in that 
category is denoted 7. 
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FIG. 26: Loss distributions for uncorrelated jump diffusion 
(solid lines), correlated diffusion with uncorrelated jumps 
(dashed lines) and correlated jump diffusion (dotted lines) 
for portfolio size K = 1 000 and 50 branches. The insets show 
the portfolio default probability PD and the unexpected loss 
UL from top to bottom for the uncorrelated case, correlated 
diffusion with jumps and correlated jump diffusion. 
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FIG. 28: The ratio between the different indicators as a 
function of the portfolio size K. The companies were divided 
into B = 1 branches and with correlation strength C = 0.5 
The indicators are calculated from the loss distribution using 
a portfolio generated by the rule set in Tab. [TTJ The ratios are 
calculated as x W /x^-V , where x is the current indicator in 
focus. 
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FIG. 27: The ratio between the different indicators as a func- 
tion of the portfolio size K. The indicators are calculated 
from the loss distribution using a portfolio generated by the 
rule set in Tab. |ll| . The ratios are calculated as a;' K - ) /x' K_1 - ) , 
where x is the current indicator in focus. 



a branch is not equivalent to a category. The concept of 
categories is developed as a tool to generate portfolios. 
However, companies which belong to the same branch are 
generally from different categories. 

The companies are assigned to the branches randomly 
for each iteration. This assures that, for example, two 
of the high risk companies in the portfolio will not be 
correlated with each other every time the simulation is 
run. In Fig. 28 the indicators are plotted as a function 



of the portfolio size K. For a portfolio of size K = 50 
the EL is improved by about 18%, and the UL by 22%. 
For a large portfolio with K = 1000 companies one finds 
the improvement in EL to be 0.74% and about 0.72% 
for the UL. The skewness and kurtosis decrease slightly, 
and the ao.999 -c fuantile and the economical capital EC 
are both improved. 



2. Drill down risk in the presence of jumps 

Finally, we examine how the jump term J(t) alters the 
pruning of a portfolio. To make this effect as clear as 
possible the correlations between the bonds are turned 
off. The jump probability is set to A = 0.01, along with 
the parameters \ij = —0.4 and a.j — 0.3. 

In this simulation, all bonds have the same parameters 
for the jump term. The portfolio is again generated by 
the rule set in Tab. [IT] In Figure |29]the various indicators 
are shown as a function of the portfolio size K. 

The benefit of removing the worst company from the 
portfolio is quite impressive. In these settings the EL is 
reduced by 10%, when having a portfolio with original 
size K = 50. The UL is improved by 9% for the same 
portfolio size. Even though these improvements are large, 
they are still somewhat lower than in the case without 
jumps. The jump term tends to smear out the differences 
between the companies in the portfolio, making them all 
behave worse. This means that the company which is 
removed doesn't stand out as much as it otherwise would. 
This fact is reflected in the ratio between the indicators 
seen in Fig. [29] 
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FIG. 29: The ratio between the different indicators as a func- 
tion of the portfolio size K, when the jump term J(t) is in- 
cluded. The parameters for the jump term are oj — 0.3, 
fij = —0.4 and A = 0.01. The portfolio was generated by the 
rule set in Tab. |ll) The ratios are calculated as a;' K - ) /x' K_1 - ) , 
where x is the current indicator in focus. 

As in the previous studies the skewness and kurtosis 
will increase when removing one company, whereas the 
ao.999-quantile and the EC will improve by as much as 
four percent for portfolios of initial size K = 50. 

VI. SUMMARY AND CONCLUSIONS 

The microscopic and dynamical character make struc- 
tural credit risk models well suited for a variety of differ- 
ent applications in physics and complex systems. We for- 
mulated a structural model for credit risk which includes 
jumps in the stochastic processes for the asset prices and 
correlations between them. As a first step, we solved a 
simplified version of the model, i.e. without jumps and 
correlations, analytically and obtained an asymptotic ap- 
proximation. Thereafter, we numerically evaluated loss 
distributions and the corresponding expected loss EL and 
the unexpected loss UL. To this end, we performed de- 
tailed Monte Carlo simulations for the full model. For 
theoretical and practical purposes, the asymmetric and 
leptokurtic character of the loss distribution is highly im- 
portant. Hence, we carefully investigated how the shape 
of the loss distribution depends on the various model pa- 
rameters. First, we illustrated that a large portfolio is 
less risky than a small one, in accordance to a basic con- 
cept of investment theory known as diversification. Sec- 
ond, we demonstrated the influence of growing maturity 
and the interplay between drift and volatility. 

Having presented the basic features of our model, we 
turned to more advanced issues which are not only of 
theoretical interest, but might also be of direct relevance 



for practitioners. The portion of the company's assets 
which is financed by loans, the leverage, affects the de- 
fault probability and the losses rather dramatically. A 
high leverage induces a high default probability and large 
losses. Since the default probability strongly depends on 
the leverage, the sensitivity of the loss distribution to the 
leverage considerably grows with the portfolio size. 

Taking the jumps into account, we demonstrated, first, 
that an increase of the jump intensity, the jump size 
mean as well as the jump size standard deviation induced 
larger losses and higher default probabilities. Second, we 
showed that the impact of the jump on the loss distri- 
bution does not depend on the portfolio size in the same 
way as, for example, on the drift and the volatility. 

Correlations heavily influence the loss distribution and 
the default probability. Although the expected loss is in- 
dependent of the branch correlation, the correlations do 
affect the features of the loss distribution. We carefully 
illustrated how the correlations hamper a convergence to 
a symmetric Gaussian-type-of shape and showed that 
the kurtosis excess develops a maximum. Moreover, the 
impact of the correlation structure depends significantly 
on the size of the correlations. Strong correlations make a 
branch act like a single company. We also showed that a 
small, highly correlated branch gives essentially the same 
model outcome as a large branch with a modest correla- 
tion. 

There is a subtle interplay between jumps and correla- 
tions. The uncorrelated jumps can partly neutralize the 
correlations. Economically, it is also meaningful to in- 
clude correlations between the jumps themselves. We 
have demonstrated that correlated jumps have a pro- 
nounced effect on the loss distribution, leading to ex- 
tremly fat tails. 

We investigated how the moments of the loss distribu- 
tion will be influenced when adding a new company to 
a given portfolio. It was shown that it is possible to ex- 
press the moments of the loss distribution for a portfolio 
that consists of K companies in terms of the moments of 
the K—l sized portfolio. In addition, we studied numer- 
ically how different indicators, such as the expected loss, 
unexpected loss, skewness, kurtosis, economical capital 
and the ao.999-quantilc, behave when removing one com- 
pany from a given portfolio. We examined three portfolio 
setups, one with and one without correlations, and the 
third one including jumps. For all three setups it was 
found that both the EL and UL was decreased whereas 
the skewness and kurtosis was increased. Furthermore, 
the ao.99g-quantile and the economical capital decreased. 
These effects were more pronounced for smaller portfo- 
lios. 

When comparing the ratios of the different indicators 
used, it was found that, in general, correlations increase 
the indicator ratios, whereas jumps tend to decrease the 
indicator ratios. 



22 



Acknowledgments 



Up to third order we have 



We thank P. Jochumzen and A. Muller-Groeling for 
fruitful discussions and helpful comments. We ac- 
knowledge financial support from Det Svenska Veten- 
skapsradet. RS is grateful for support from Deutsche 
Forschungsgemeinschaft under grant no. SCHA 1462/1- 
1. 



APPENDIX A: LOSS DISTRIBUTION FOR A 
HOMOGENOUS PORTFOLIO 



We express the distribution (17 1 in terms of its Fourier 
transform, 



+ OC 
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where (L k ) is the n-th moment (13 1 of the distribution 
p k {L k ). This yields 

Khx{{l-P D ) + P D Q{u>/K)) 
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(Al) Collecting everything, we arrive at the expression ( 19 1 



where 



"M = II / dL k p k (L k ) 
k=i{ 

J dI k p k (I k )exp(iuiL k I k /K) . (A2) 



is referred to as characteristic function. If we assume 
that the face values and the parameters of the geomet- 
ric Brownian motion are the same for all companies k, 
the distributions p k (L k ) for the individual losses become 
the same for all companies. This also holds then for the 
distributions p k {I k ) of the default indicators. Hence, we 
find 



K - 

r ( W ) = II / dL kPk(L k ) 

k=l{ 

((1 - P D ) + P D cxp {iujL k /K)) 
= (J dL kPk (L k )((l-P D )+ 



K 



P D exp {iujL k /K)) 

exp(inn ((1 - P D ) 

+P D Q(u/K))) 



with 



Q(lu/K)= J dL k p k (L k )exp(iujL k /K) 
o 



(A3) 



(A4) 



As we aim at an approximation for large K, we expand 
the exponential to obtain an asymptotic series in 1 / K. 



APPENDIX B: LOSS DISTRIBUTION FOR AN 
INHOMOGENOUS PORTFOLIO 



We generalize approximation (19) to the case of an 



inhomogenous portfolio with varying face values F k and 
individual loss distributions p k (L k ). To this end we write 
the portfolio loss as 

i=%^ = E^, (Bl) 
Z^fc=i ^fe fe=i 

where ^ k = Pk/^2f=iFj is the fraction of money in- 
vested in obligor k. Following App.[A] the characteristic 
function is then 



r(w) = p (i 1 - p D,k) + Pn,kQk{u)J , 

with 

l 

Q k (ui) = J dL k p k (L k ) exp (iuj-f k L k ) . 
o 

Expanding the exponential in Q k (u>) yields 

Q fc (w) = 1 + iu; lk (L k ) - -uj 2 j 2 k {L 2 k ) 



(B2) 



(B3) 



6 



(B4) 



which finally leads to 



P{L) 



1 

27 



dui exp — iu> 



CX P i-^illl(PoALl)-P 2 DM {L k f) 



k=l 



23 



K 



fc=l 



+3Pj> {L k ){L\) + 2P% JL h 



(B5) 



/+oo r+oo 
dhpi(h) ■ ■ ■ <IIkPk{Ik) 
-oo J — oo 

fF 1 rF K 

x / dTipi(ri).-- / dT KPK {T K ). 
Jo Jo 



as a generalized approximation for the loss distribution. We introduce the sum over all the face values 



APPENDIX C: DRILL DOWN RISK 

We want to express the moments of the loss distri- 
bution for a portfolio which consists of K companies in 
terms of the moments of a portfolio which consists of 
K — 1 companies. The loss of an individual company 
is given by Eq. ( 25 1 and the total loss for a portfolio 



K 



containing K companies by Eq. ( 26 1 . To obtain the dis 



tribution of the entire portfolio loss, one needs to average 
over all distributions of the individual losses pk(Lk), and 
over the indicator distribution Pk(Ik)- This is done in the 
following way, c.f. Eq. (17 1: 



dhpi(h) 



+oo 



dTiPi(Ti) 



dI K PK(jK) 

dT K p K {T K ) 

(CI) 



where the upper index (K) indicates that the total num- 
ber of companies that are averaged over is K. 

The moments of the loss distribution arc defined as 

((l {K) Y) k - J 1 (L^y i K \L^)dL^ K \ (C2) 

where L^ K ^ is the loss for a portfolio containing K com- 
panies, p( K \LS K }) is the corresponding probability den- 
sity function and n is the order of the moment. The sub 
index K in the bracket indicates that the distribu tion 
for K companies, p^ K \ is used. Plugging equation (CI) 



into (C2) we find 



K»)"> 



K 



(C3) 



dhpi(h) 



-oo 



+ OG 



dlKPK^jK) 



dTiPi(ri)- 



oc- 
Fk 



dT K p K (T K ) 



x 8 iW 



Integration over dL^ yields 



(C4) 



3=1 



(C5) 



and rewrite the sum in Eq. ( C4 1 as 



a-i) 



try E r ^ • ( C6 ) 



3=1 



Using the binomial theorem to rewrite the n:th power 
sum, we find 



A' 



K-l) 



< + 00 



/-t-oo /-t-oo 
dimih)--- / di K p K {i K ) 
-oo J —oo 

[ 1 dTipi(ri)--- / Va-paOa) 

Jo Jo 

(\ n— f 



fc*"- 1 ) 



(C7) 



This can now be expressed in terms of the moments of 
the corresponding K — 1 distribution, ^(Z^ -1 ))™^ , 

and the moments of the new distributions, pk(Ik) and 
Pif(Tif), in the following way 



K>)">, 



S (;)( (£ <-)' 

(^(A-l))*, ' 



(C8) 



where i/ is an integer. It is important to remember which 
distribution the bracket notation (• • •) refers to. The mo- 
ments of Ik and T k are taken over the distributions 
Pk{Ik) an d Pk(^k) respectively, and the moments of 
Ij( k -i) are tajjgn over the distribution pi K ~ 1 ^ (ZA-^ -1 )) . 
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